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ABSTRACT 


A new, detailed first principle kinetic theory for electrons is 
presented which is neither a classical fluid treatment nor an exospheric 
calculation. This new theory illustrates the global and local properties 
of the solar wind expansion that shape the observed features of the 
electron distribution function, f , such as its bifurcation, its skewness 
and the "differential" temperatures of the thermal and supratherraal 
subpopulations. Our approach starts with the Boltzmann equation and 
retains the effects of Coulomb collisions via Krook collision operat 
without recourse to wave-particle effects. We conclude that 
Coulomb collisions determine the population and shape of f of both the 
thermal electrons energies (E < kT) and suprathermals (E > kT)., We find 
that electrons with E > 7 kT constitute a special subpopulation of the 
suprathermals — insofar as Coulomb collisions are concemed—these we 
call "extrathermals", The electrons in the thermal portion of f have 
undergone ~ 10-20 Coulomb collisions for cumulative momentum transfer 
en route to the observer at 1 AU; this population is thus more removed 
from the properties of coronal electrons than is the suprathermal 
population. This latter group retains a strong memory of coronal con- 
ditions, since they have undergone only a few momentum transfer collisions. 

The thermal population is most nearly in collisional contact with the 
local dynamics (compressions, rarefactions, etc.) of the solar wind. The 
suprathermal portion of f is determined by Coulomb collisional inter- 
actions with the distribution of solar wind material on radial scale of 
the Heliopause itself. In this respect the suprathermal portion of f g is 
found to be responsive to the consequences of the global dynamics of the 
solar wind expansion. We find that this subpopulation is an attenuated 
vestige of collisional populations deep in the corona, (1.03-10 1^) , which has 
been redistributed via Coulomb multiple pitch angle scattering on magnetically 
open field lines. The suprathermal particles moving toward the sun are computed 



to be observed as a result of Coulomb collision induced backscattering at 
larger (1-10 AU) heliocentric distances than that of the observer. Based 
on this theoretical picture quantitative estimates for the partition of 
thermal and suprathermal phase density, the break in the velocity dis- 
tribution and the magnitude of the skewness (heat flux density) agree 
well with those typically observed near 1. AU. These calculations predict 
that the extrathennal fraction of the phase density, the extratlieraal 
temperature, and the net heat flux density carried by electrons should 
be anti-correlated with the local bulk speed in quasi-steady state flows, 
and that the radial variation of extrathermal temperature 

inside 1 AU should be essentially independent of heliocentric distance. 

Our work also shows that the observation of suprathermal particles 
cannot be taken as a priori evidence for in situ wave particle interaction 
since we can theoretically calculate a suprathermal population of solar 
wind electrons at 1 AU by assuming wave particle interactions are not 
present any where in the heliospheric cavity. This new way of obtaining 
suprathermal particles arises from the consideration of the non-vanishing 
polarization electric field that is always present in an inhomogeneous plasma 

This work suggests that the local cause and effect precept which 
permeates the physics of denser media must be relaxed for electrons in 
sparse and radically inhomogeneous plasmas such as are found in the solar 
wind between the lower corona and the interstellar medium* The local form 
of transport laws and equations of state (e.g,, Q = -HVT, P = NkT) , which 
are familiar from collision dominated plasmas, must be replaced with 
global relations that explicitly depend on the relative position of the 
observer to the boundaries of the system. 
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I NTRODUCTION 

In this paper (I) we suggest that the typical distribution function 
of solar wind electrons in a quasi-steady flow is shaped primarily by the 
properties of Coulomb collisions and the "smooth" macroscopic fields 
implied in the solar wind expansion. Following a straightforward 
mathematical procedure, we illustrate how some of the controlling effects, 
both local and global, are reflected in the in situ observations. 

In companion papers we discuss: (II) the consequences of this theory 

and its comparison with known experimental facts, (III) the impact of 
this approach on the question of transport phenomena (such as heat flows 
viscous effects, etc.) in the heliosphere and the relation of the current 
formulation to canonical transport methods, and (IV) a mathematically more 
rigorous formulation of the problem discussed in I. 

Historically, theoretical treatments for solar wind electrons have 
variously approximated the electrons as either: [1] part of one fluid 

model of the wind (Parker, 1963), [2] as a separate fluid interacting 
via Coulomb (local) momentum transfer with the ions (Hartle and 
Sturrock (1968)), or [3] as an exospheric ("collisionless") population 
above some radial distance (Jockers, 1970)). Ogilvie and Scudder (1978) 
showed with in situ measurements between Q.45 and 0.9 AU (Feldman et al . 

(1978) at 1 AU concur) that the local mean free path, X, for scattering of 
thermal electrons (core) via Coulomb effects was less than the conventionally 
defined scale height, H, implying that Coulomb effects could not be 
neglected even within 1 AU. Therefore, the propagation of electrons 
neglecting particle-particle collisions from a baropause to the in situ 
observer is not the correct kinetic description of the solar wind’ s 
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expansion insofar as electrons are concerned. However, Ogilvie and 
Scudder (1978) indicated that X/ll was not very small compared to unity 
indicating that the fluid treatment was not justified either. It is 
therefore fruitful to explore the effects of Coulomb collisions in 
shaping the in situ distributions, by not making either the exospheric 
or fluid approximation. 

The observations of Ogilvie and Scudder (1978) (Feldman et al. (1978) 
also indicated that the local scattering mean free path (mfp) for "halo” 
electrons was long compared to the local density scale height. In the 
customary sense of the word, the halo particles were locally "collisionless". 
Nevertheless, there Is some weak, omni-present, unavoidable Coulomb Interaction 
with the thermal plasma. The suprathermal electron population of fixed energy 
is generally observed in all octants of velocity space with comparable 
phase density Montgomery et al., 1968, Scudder, 1970, Ogilvie et al., 1971. 
Local wave particle scattering has been invoked since the discovery of the 
solar wind to explain the near Isotropy of all plasma components in their 
proper frame. Feldman et al. (1975), for example, have suggested this as 
the reason for the "ubiquity" of suprathermal electrons. 

The origin of the spectrally distinguishable halo subpopulation of 
solar wind electrons was Initially suggested by Ogilvie and Scudder (1978) 
when they extrapolated the electron measurements made by the Mariner 10 
Plasma Science Experiment of the radial variations (0.45 < R < 0.9 All) 
of the core and halo (differential) "temperatures" 

T = 


L f dlnfyi 
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to determine the radial positions where these best fit variations 
predicted a common temperature. Depending on the assumptions this zone 
of common (differential) temperature was between 2-15 1^. The suggestion 
of that paper was that a fractionation process of the antecedents of halo 
and core subpopulations had begun in this radial range. The subsequent, 
approximate experimental inference of the origin of the u strahl" portion 
of the halo in very high speed flows (Feldman et al . , 1978) confirms the 
approximate location of this fractionation process. 

We shall demonstrate below, by explicit calculation, that the supra- 

thermal electrons with E > 7 kT observed near 1 AU are an attenuated 

c 

vestige of collisional populations deep in the corona, (1.03 ~ 15 R^) , 
which has been redistributed via Coulomb multiple small angle pitch angle 
scattering. The cumulative, spatially integrated effects of Coulomb pitch 
angle scattering are sufficient to explain the observed proportion of supra- 
thermal electrons both fore and aft along the local magnetic field. We 
show that the halo population is shaped by the spatially integrated, 
non-local (x.e,, globa]) features of the heliospheric cavity — such as its 
size, its large scale magnetic topology and the spatial variation of its 
contents — especially near the extremities of the solar wind expansion. 

In this latter sense in situ electron observations allow remote observations 
of important regimes of the solar wind expansion which are presently not 
accessible. 

Terminology 

It is appropriate to comment on our use of the word "suprathermal" 
and other new terms we wish to use in the remainder. From the lntin 
su P ra > meaning above, the clear intent of 1 'sup rathe imal" is to delineate 

(where k_ is Boltzmann 


those particles with energy larger than E » k^T 
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constant which is the kinetic energy of a particle moving with the most 
probable speed. For electrons* we operationally assign the ''core" 
temperature to delineate the most probable speed. For particles with 

E> k-T it is consonant with the literature to discuss this entire 

T5 core 

range as the suprathermal domain. In the companion paper, II, we shall 

O-r 

show that electrons with kinetic energy E > 7 kT c are a special subset 
of the suprathermals insofar as Coulomb collisions are concerned. These 
particles have energies sufficiently beyond that of the thermal electrons 
as to have a minimal local Interaction with them. Moreover, as shown 
below these extrathermals have also come from substantially beyond the 
I sition of the observer. For clarity we have labeled this group as the 
"extrathermals" noting that extra means "beyond” in Latin. The particles 
of intermediate energies (1 < E/k^T^ < 7) are most suitably labeled as 
"transthermals" (trans = across) for two reasons: Cl) they bridge the 
energy interval between sub-thermal and extrathermal populations and 
(2) these electrons are populated from collisional antecedents with 
non-negligible contribution distributed across a substantial portion of 
the heliosphere. We shall use the terms subthermal, thermal , trans- 
thermal, and extrathermal to delineate specific energy regimes in the 
above sense? when we do not wish to distinguish between trans* and extra- 
thermal we will Use the word suprathermal. As discussed below we do not 
regard the term, "suprathermal" as synonymous with that suggested by Parker 
and Tidman (195S) for particle populations generated by local Fermi 
mechanism and/or betatron action. 

Organization 

We outline the structure of this paper to orient the reader. The 
theoretical transport problem is posed in a formulation section, Cl). 
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We next discuss formal solutions, (2), which give the flavor of the explicit 
solutions whose details can only be obtained after the mathematical con- 
siderations, (3), and implementation section, (4). The outcome of two 
"realizations” of the formal solution are given in a results section, (5), 
followed by a discussion, prediction and conclusion section, (6). 
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1. 70KHUIATI0N 

For the present we desire to estimate in the simplest possible 
fashion the consequences of a physically complete picture of Coulomb 
collisions as they would shape and determine in situ observations of 
electrons in the solar wind without considering wave-particle scattering. 

The advantages of the present formulation include simplicity and ease of <- 

physical interpretation. Estimates of the importance of Coulomb inter- 
actions can be made simply by evaluating integrals * rather than solving 
partial differential equations. In this approach the solution of the 
scattering problem is available in closed form and has a conceptually 
simple interpretation with predictions which can be checked against the 
variability and correlations of local observations! it also - allows definite 
predictions at other heliocentric distances. 

The basic equation governing our discussion is the Boltzmann equation 


7 f £ 2 f 

at ^ V' d* + 


Sf 

Si 


^ Collision 


CD 


where the symbols have the usual meaning; f is the particle distribution 

function at coordinates (r» v) in the phase space; _F are the forces ex- 

6f 

perienced by the particle and is the Boltzmann collision term which 

is given explicitly by 


if I 

'Collision 


n [ fj&fdJL isi 4 tv) ( v.y, j-f ry, ) 

- Jdv'fda' 4(y r ) if-v, 7/cy/)], 


wh,ere'V_, are the scattered velocities of these two particles* 


The first term of 


determines the loss rate of f from 


collision 

(r., v) as a result of the scattering of particles with pre-collision 


ORIGINAL pagbjS 
of POOH QUALITY 
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velocities v’ , v^' . The differential scattering cross section is given 
hy da /dQ. 

For purposes of the estimates in this paper we shall approximate 
the collision term in the Boltzmann equation in the form proposed by 
Gross and Krook (1956) ; this approach is also known as the relaxation 
time approximation (c.f., e.g. , Reif, 1965); the approximation is 
mathematically exceptionally simple and physically not only instructive 
but has been shown to be quite satisfactory in other comparable physical 
situations. 

As clearly realized by Gross and Krook (1956) , this approximation 
is not restricted to the small mean free path transport case which leads 
to the Chapman-Enskog type analysis. Therefore, our use of the relaxation 
time approximation is not synonymous with the short mean free path limit 
(vt/L << 1 , where v is the particle’s speed and L is a characteristic 
scale length of the gradients). 

Consider at some point jr, the (unknown) electron phase density 
f(£, R, t). Then in Krook’ s approximation, f obeys the equation 


At 


v VI 

'Z'C'c, v) 


( 2 ) 


with the usual definition 



and where E_ and are the macroscopic ’’smooth” electric and magnetic 
fields and jg is the acceleration of gravity. (Remember we assume no 
wave-particle interactions.) By comparing the structure of (2) with 
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that of (1) we obtain the correspondences: 

ZT = ( nfdv, fdJl i?S(y.) Iv-v.l) 

j: . ^(JdVi'JjJL' ah' ■f(y')-fcy'i') i v'-y ! ) 

- — — ■ — 

(]<**, Jssi in fci,)i?-y.i) 

where it is clear that and 7 are functionals of the unknown distribution 
function. 

The Krook approximation to the collision operator 6f/6t not only 
retains the formal structure of the Boltzmann operator, but allows, by 
a judicious choice of (r:, v) , a linearization of the full transport 


equat ion 

- 6f/ , 

of 1 6tJ 


Krook 


instructive. We shall now examine some properties 
to understand the role of f. (r , v ) in this approximation. 


In the special circumstance of a uniform system where electric 
fields and gravity are unimportant and the distribution function is 
perturbed isotropically in the proper frame velocity space and uniformly 
in configuration space, then (2) reduces to 

4? = = 

dt di. rg 


(4) 


The solution of equation (4) is 

JYwtW- + Q(m) -exp(- (5) 

where for this example t is assumed independent of time t. From this 
example it is clear that should be interpreted as a final state of 
the relaxation process. The temporal relaxation given by (5) is caused 
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(a) by (f -£.)! r 0, (i.e., the departures of the initial from the 
“■ to 

final state) and (b) by the evidence for a finite relaxation time, T- 

The preceding discussion was for an extremely special situation for 
the Boltzmann evolution operator with a collision term approximated 
a la Krook. By contrast we require the steady state solution of Equation 
(2) = 0) for the more general situation when there are non- 

negligible forces and spatial inhomogeneities in the expanding solar 
corona. This circumstance is far more common than the simplifying 
example of (5). Nevertheless that example has been instructive in 

supplying insight into the meaning of f which may be carried over into 

K. 

the more general case we wish to pursue: namely if f (r_, v) = (r_, y) 

then this is a stationary sroiution of the Boltzmann equation with a Krook 
approximated collision operator. 

The distribution f^(r, v) represents the time independent, spatially 
inhomogeneous final state towards which the relaxation mechanism, with 
scale time T, drives an initially perturbed distribution. Of course, f, 
is not known a priori ; however, using the in situ tendency for the thermal 
electrons to be nearly Maxwellian (Montgomery, 1968, Feldman, 1975, 

Ogilvie and Scudder, 1978) we make an empirical ansatz for the final 


profile in seroth-order approximation: 
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fKl e \V~IAtr 


which is a connected, Maxwellian distribution with spatially varying 
density, n(r), hulk velocity U(r) and temperature, T(r) . 
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We note explicitly that this estimate (suggested by data) is the 
initial step in a chain of iterations: 


d£°> 

<±6 


M. 


^cfn 


(7a) 


where 


d £ 

d-fc 


M) 


TCCl 


m 


(7b) 


t* = t 


( tt -0 


( 8 ) 


where we have denoted functional dependence by the square brackets. In 

this paper we shall consider the solution of (7a) and we shall discuss 

the impact of subsequent iterations 

The relaxation process visualized here is characterised by the 

familiar rate of Coulomb collisions v between "test" electrons and 

c 

"targets", i.e. the positive Ipns and other ambient electrons. Recalling 

that v s 1 /t we have (c.f. e.g., Rossi and Olbert, 1970) 
c 


i 1 ^w T J e J( w kO 

K " Jw? ' 


(9) 


This is the rate at which phase density is scattered out of the observer 5 s 
phase density element, where in cgs units 

, n ^CiW(T''*T e yn''0 

' 0.3S ‘ 
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and where T is the harmonic temperature of T g and the speed jwj = w 
is the particle's speed in the proper frame. The root mean square (rms) 
thermal speed of species a is given hy 




(k B 


Boltzmann 1 s constant) 


The statisical effects of the thermal spread of the ambient targets is 
contained in the function r, of (9) defined by: 
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A more general form for-- is given in Rossi and Olberc (1970). In 
that expression *' . is a functional of the ambient electron and ion 

distribution functions f^, and f^. The expression given in (10) is 
evaluated consonant with 7(a) as vC (f- f, - ^) s where we have 

"" S K j 1 

used local Maxwell ians for both electrons and ions. 
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2. INFORMATION FROM FORMAL SOLUTIONS 

Although formal solutions are rarely directly amenable to calculation, 
they often compactly display the conceptual character of the more com- 
plicated explicit solutions. We proceed to give several instructive, 
equivalent, formal solutions to the Boltzmann equation, with the Krook 
approximated collision term, in order that the reader obtain the flavor 
of the more complicated explicit computational form of the solutions which 
we discuss in the next two sections. 

The formal solution of (2) at the observer’s time t is 

J expc-s (*.*'))£ If l/nr] di' a 


where 



( 12 ) 


and where the indicated path integrals ^ are performed following the 
trajectory of a representative point of an electron in the six-dimensional 
phase-space. 

Noting that S(t, t) = 0, and that 


-=• 


-_dt' 


V ~t'<i 


and, with the definitions 


(13) 


C14) 




and 45 = Sffjt'J 
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equations (11) may be recast in a number of illuminating ways: 


fit) = f^ltf) exp(So) *J? exp(- /)s) 0\S)]) d(fc) (15) 


which, upon integration by parts, becomes 


lit)* S t (&S=o) -tj expOtS)^ ~3 q&) 


A r -t'w] ) . 

y das) 


J 


(16) 


or, equivalently 
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d£ 

Set ) - S k (M ~ d ) + jf dp (- tsfc d£{] 


J* f6S*o) 


(.17) 


Equations 11, 15-17 are equivalent, highly compact, formal answers to 
the transport question posed at the beginning of this section. Before 
discussing the content of these solutions we need a physical interpretation 
of the quantity AS. 

The quantity AS is the collisional analogue of the optical depth 
of radiative transfer. We will refer to this intrinsically non-negative 
quantity as the collisional depth of the observer relative to the source. 

It is a cumulative Qj:) counter of collisional cycles (1 /t) experienced 
by the representative point of a particle in the phase space element 
during its motion £at t ! 'j from the departure at t’ to the present, at t. 

The quantity exp (-AS) which occurs throughout these equations can be thought 
of as a normalised survival probability upon transiting the collisional 
depth AS. 
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To appreciate some of the essential physics of the non-local terms 
of (15-17), we must discuss electron kinematics and the use of collisional 
depths and boundary conditions. 

The electrons most probable speed throughout the solar cavity is many times 
the typical solar wind speed. Hence, nearly all electrons are kinematically 

able to transport information in both directions along a magnetic tube of 
force. Because of this mobility, it is by no means obvious where the 
electron has been in the solar system preceding its detection. For 
electrons the in situ observation from a given collisional depth AS of 
the velocity distribution function is a synthesis of particles and in- 
formation coming from both fore and aft along the local magnetic field. 

Collisional depths, like their optical counterparts, are 
intrinsically non-negative. When considering contributions to the local 
distribution from a given collisional depth AS, we are considering con- 
tributions from particles which at time t’ < t were either inside or 
outside of the observer’s radial position (c.f. Figure 1). A corollary 
to these comments is that all electrons observed at _r at the time t 
probably did not leave the sun at, or even nearly at, the same time, t^, 
in the past. Rather, the observed f(r, V, t) is the superposition 
(i.e., an integral) of all the surviving probability densities, whose 
representative electrons have accessible trajectories to the current 
observer's point (r, v) in the full six dimensional phase space. Depending 
on the collisional depth and its variation with energy and pitch angle, 
and the temporal boundary conditions at the source layers, f(AS (t T )), the 
observed distribution will be more or less reflective of the coronal or 
collisional conditions along the local observer’s magnetic tube of force. 
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Thus, equation (15) has the following physical content: the 

velocity probability distribution f(r, v, t) observed at time t, is the 
sum of two parts — the first on the RHS of (15), the so called boundary 
term, is the contribution to the currently observed f that has survived 
extinction from the space time boundary of the system to the observer. 

The second contribution is the folded sum (J:) over all intermediate 
collisional depths 0 < AS < S of f,[t' (AS)], weighted by their 
respective attenuations, exp (-AS), between source and observer. 

(Remember that for a given AS there are two positions along the tube 
of force from which $ has contributions: one effective location outside, 

one inside that of the observer, c.f. Equation (35).) 

Equation (16), although entirely equivalent to (15), gives a 
slightly different physical description to the observed distribution: 
the first term on the RHS of (16) is the "local" term, since it is the 
source distribution at collisional depth zero from the observer — hence 
the name; the second term describes the corrections to the local sources 
as properly attenuated signals of the variability of the non-local 

distribution with differential changes in the collisional depths from the 
observer. This formal solution (16) is the point of departure for our 
discussion of III. Equation (17) follows from (16) by change of variables - 
Note from (17) that if the system is everywhere homogeneous, then 
f^(AS = S q ) = f^CAS = 0) and (17) collapses to f = f^(AS). This is the 
special case of thermodynamic equilibrium. Conversely, if for example, 
the electron temperatures at AS *= S^, 0 are not equal then the J: term of 
(17) is not equal to zero (since the integrand is positive definite). The 
corrections to the local population (f^CAS =0)) that are observed are the 
signatures of spatial inhomogeneities of f^ in the system. 
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From these formal solutions we see that the local observations 
represent a superposition, of phase densities displaced from the observer 
(AS^O) folded with the probability that these phase space elements can 
"survive" between their collisional origin and the observer. This is a 
mathematical statement of the dynamical accessibility of the electrons to 
the observer. In order to further evaluate these solutions some 
mathematical arguments are required. 
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3. MATHEMATICAL INTERLUDE — CAUSALITY 


The discuss ion which follows pertains to the directions of 

information flow to a given point in (r, v) space. 

In order to t, reasseInble ,, the Lagrangian style solutions available 

from. (15) for the evolution of different phase space elements into an 

Eulerian description for the velocity distribution at a given position 

r, as t : we must categorize classes of trajectories in steady state which 

are accessible to the position of the observer. 

In this connection it is important to note 'that electrons of energies 

less than 10 keV have gyro radii which are extremely small by astrophysical 

scales at any heliospheric distance. It is thus difficult for solar wind 

electrons to drift perpendicular to the magnetic lines and therefore, 

excepting the effects of collisions, the electron guiding centers reside 

on the line of force on which they were initially injected. Although 

collisions can enhance crossfield diffusion, this effect is negligible 

— L6 

at r = 1 AU since D /Dij ^ 10 where D it are the collisionally induced 

o -l 1| -Ml 

diffusion coefficients transverse and longitudinal to the local magnetic 

field. This ratio scales as ORIGINAL PAG0 IS 

OF POOR QUALITY 
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= 17 I&> ) 3k jffeogto 

Di/Uu l D L I Al re) 
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-24 , . , 

This ratio is approximately 10 in the corona and is small throughout 
the heliosphere. Thus, the history and future of the guiding center 
motion of electrons in the solar wind is intrinsically one- dimensional and 
determined by the topology of one given line of force. Nevertheless, the 
guiding center motion does not behave like collisionless theory would 




MAGNE 


Figure 1 
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indicate; rather, there are cpilisional modifications within this one- 
dimensional potion. Accordingly, we adopt an archimedian magnetic topology, 
Figure 1, after Parker and use the arc length, JL*, along the given lines 
of force interchangeably with the corresponding radial distance, r*, of 
the observer from the sun's center. As far as the motion of electrons is 
concerned, the history and future of the electron trajectory is along the 
given tube of force that passes through the observer’s position at (r*, t*) . 
The discussion which follows assumes the fluid flow to he time independent 
in the corotating frame. In paper XX we will discuss to what extent the 
violation of these assumptions will change our results. 

We introduce some diagramatic methods for listing the types of 
particle trajectories that can gain access along a single tube of force 
to the observer 1 s phase-space volume element, drdv about (_r*, v*) . In 
Figure 2 we illustrate representative trajectories of four types of con- 
tributions to the locally observed distribution function, f. There are 
two classes of contributions to the forward pitch angle portion, f , of 
f: (1) those contributions with ultimate ’’source” in a collisional 

distribution outside of the position of the observer, (&*), for times 

(t^) prior to the time (t*) of detection; and (2) those contributions from 
inside (& the position of the observer at times t^ prior to the detection 
time t*. (Refer to Figure 1 for the location of and } These same 
two categories apply to the available contributions to the aft, f , pitch 
angle distribution. 

The first group of accessible trajectories in Figure 2 are those 
which leave the source point (1) with space-time coordinates (& <5 tM 
and arrive at the observer’s point labeled (0) with coordinates (&*, t*) 
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without having been pitch angle backscattered. We introduce a special 
symbol for this type of collisional depth between source and observer, 

6S^q* Note, however, that in general OS f- AS. If we desire to consider 
contributions from all collisional depths, even from this restricted 
type of contribution from this one source layer, (t* -t^.) must be allowed 
to range from 0 to « in order to properly compute these contributions. 

By definition then the probability, P, of transiting is given by = 

exp = P° exp (-6S^p), which defines s 1. 

Another example of this group are those trajectories which are 
like 2 1 0 In Figure 2, characterized by having one stochastically 

induced turning point between source layer and observer. Collisional 
depths are additive (c. f. , (12)) and therefore 

&S zo - flS a( -J fl.S l0 = AS 2l +^S /e ag) 

Therefore , = exp (- /\S zl ) <?Xp C- $S, 6 ) ^ T^eicpC- <£$*« ). 


Clearly, AS^g may also range from 0 - c 3 , either because AS^, 
AS 1Q do. Similar considerations, yield 


and/or 




( 20 ) 


with any and/or all the quantities having infinite range. Higher order 
quantities of these types can readily be seen to be more turning points 
between the time and particle left the source and reached the observer. 

We shall assume for paper (I) that elastic pitch angle scattering 
is more efficient than Inelastic scattering in backscattering electrons 
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everywhere in the heliosphere. In the present context elastic scattering 
implies that the test electron undergoes a pitch angle change with no 
loss of proper frame kinetic energy; inelastic scattering implies a loss 
of proper frame kinetic energy while the pitch angle is changing. 

The dominant fraction of momentum transfer collisions for test 
electrons is off of ambient ions. This is an essentially elastic process. 
Test electrons scatter off of other background electron distributions both 
elastically and inelastically. 

This elastic scattering approximation is most accurate for extra- 
thermal electrons and will be relaxed in Paper IV. The principle impact 
of this approximation is to limit the precision of the predictions of 
the sunward propagating extrathermal phase density in the trans and 
extra thermal regime. The predictions of the current formulation for 
the outward propagating extrathermal population are essentially independent 
of this assumption. 

As will be shown below, if only elastic scattering is considered 
then as shown schematically in Figure 3a 

/• / i 

JMh ___ — ~ 

to). -9 « -J- +(i m ) 

independent of energy in the extreme extrathermal regime. The super- 
scripts on f denote cuts of distribution function along (+) and, opposed 
C-) to the outward ray along the tubes of force. In support 

of our theoretical arguments that elastic collisions should be more 

important than Inelastic in determining f ( h raw ) , we discuss an ex- 

+ “ 

perimental determination of f /f as a function of energy in Figure 4. 

Here we have plotted the ratio. 
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* 


versus parallel proper frame speed, jWjjj for a published example of 
IMP -8 data (open circles) . We have indicated the theoretical asymptote 
of which follows from ignoring Inelastic collision. The experimental 

data for this published example does asymptote in the vicinity of the 

1 

theoretical value of 7 . (The solid circles show the trends of a 

6 

'’realization" discussed in section (5).) 

If inelastic Coulomb scattering assists extrathermal backs cattering 
via the speed dependence of the cross section, then 


/7 ;m<) J + 


(ZK,C ^ 
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This Implies that the differential slope of the extrathermal population 
going towards the sun will have cooler differential temperatures than 
those going away, (Figure 3(b), (c)), if Coulomb collisional losses are 
a substantial factor in the extrathermal backs cattering history. In 
the extreme limit of completely inelastic backscattering the extra- 
thermal temperature of f may even be cooler than that of the local 
thermal population (Figure 3(c)). 

The exclusive elastic case of Figure 3(a) is the most conducive 
to fitting the halo subpopulation with a drifting bi-Maxwell ian, where by 
construction 1 


dE _ AJf 




as is implicit in the modeling done by Feldman et al. (1975) for a wide 

2 

variety of data. Nevertheless the average normalized X for these 
parameterizations is rather high indicating a probable admixture of 
some inelastic backscattering (Figure 3(b)). In this latter circumstance 



24 


the least squares model of a drifting bi-Maxwellian is no longer 

sufficiently flexible to accommodate this spectral shape. If a 

bi-Maxwellian is forced to fit data as schematised in Figure 3(b) it 

2 

will result in a larger normalized y than data as schematized in Figure 
Figure 3 (a) . 

With the above provisos we shall proceed on the assumption of 
exclusively elastic Coulomb backscattering (Figure 3(a)). By simple 
geometrical considerations and the assumed stationarity of the 
Hamiltonian, we obtain 


and where 6S ^ is the collisional depth between the m’ th reference 

point and the next turning point on the trajectory as time increases. 

Therefore, the integrated probability that a particle leaves 

(& , t ) and returns to (i = & , t ) , along the tube of force of 
m m n m n 

infinite length, after only one turning point, is 
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( 22 ) 


using C21) we obtain 

P (t) - cliffy) exp(- $,x)( * ~ = Z 


( 23 ) 
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( 2 ) 

Similarly the integrated probability P that a particle leaves 
and returns to the same level after two turning points is 


V 


< 2 ) _ pCO 


• P C,) * ( 


i) 
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since the cumulative survival is a joint probability. By induction, the 
probability of surviving n turning point trajectories between a source 
position Jl* and the same source position at a later time is 

[if « 5 ) 

Thus the contributions to exp (-AS) (AS) on the RHS of (15 ) from 
an "emission” layer at for all times t prior to t* which is at 
collisions! depth SS to the observer is 


p j + o.(eU*) (26) 

=■ | (ss i-e'j) ■ <*pc- a ) + « a;> £*) 

where Q(&^, Z*) are the contributions from collisional distributions at 
Z> larger than Z*. In deriving (26) we have made use of the fact that the 
forces considered were conservative, that f^ is isotropic, and does not 
depend explicitly on time. 

In Figure 2a we also show some of the leading contributions (4- ~+ 0): 

(5 -* 4 -+ 0) to the additional term (Q) in (26) which are those contri- 
butions to the observed f from the same optical depth AS but source 
regions beyond that of the observer. In order for these source regions 
to have access to the same phase space element dr*dv* about (r*, v*) as 
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the contributions explicitly given in (26) , at least one turning point 
en route from a collisional distribution to the observer is required. 

Thus, the probability sum for this additional term is exactly that already 
used minus the first, or direct ascent (descent), term (P^°^ = 1); there- 
fore, the complete expression for (26) for all contributions from 
collisional depth AS, for particles with a specific positive d&/dt = v, ( , 
and energy E at (&*, t*) is 

eq>c-as)£ ( .t£Y5)j = <?*j> C-<Ts)-^ (27) 

For particles at moving toward the sun (i.e., negative cosines 
of their pitch angle) (Figure 2(b)) direct access to the phase space 
element (r*, v*) is not possible from inside; rather, only from outside; 
reversing the weighting of (27) we obtain 


expc-(\s)l t (as)j =. j 

When (27) and (28) are normalized by the condition 

• t 

Aa expC- j «. 4 k 0&-Q) 

&S-X) 


(28) 


(29) 


the equations (27) and (28) become 


€xp(-te)^te)!'= e<p£-fc)J| £(*,[«]) + 


(30) 
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and 


■xp£- - exp C-SS)1 j f k (t C*S]) * | i (*>CtSl)j (31) 


Equation (15) then becomes explicitly 


) r/i 


r /t\ 


\ * 
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(32) 


; ’ i f eK F c-K[t^i) f k c±'(e'(rn )]tt Jj'jy 
\ ?/ {/*>) r &(*'&))] dt dr 1 

*\ r 

+ ' fe ! / eXpC-JSU-t'O-' }])+!< [t V^Vr^)] ' di* df' 

v 1 i \th) T[t '(*'(<■<)}] dr dr' 

From the Parker topology of the interplanetary field defined by 
bulk speed U, solar rotation rate 0, at location r T and ecliptic colatitude 
, 8 1 , tore have 


dr 1 


V 6¥ + u 




(33) 


We note that r(£*) £ r (t ) and that 

u o 


£ - + 


Uj2' 


dt* 


d£‘ 

dfc' 


in the first integral of (32) 


and that 


cR} 1 

d-i 1 


cR 1 


in the second integral of (32) 


By conservation of energy and magnetic moment 


{ V n 2 (r‘) = V t( L th) + 





where 0(r) is the interplanetary electrostatic potential and m V M and 

e II 

m e V ± are the particle’s components of momentum along and transverse to B. 
The gravitational terms in the conservation of energy have been neglected. 


Thus (32) may be computed as integrals over ordinary space to be: 



' n (*({«]) 


The most prominent advantage of (35) is that all explicit references to 
turning point orbits are absorbed into the coefficients of 1/7 and 6/7. 
The remaining quadratures are for unscattered orbits between collisional 
depth 6S and the observer. It is this circumstance which allows the 
use of collisionless guiding center relationships, (34), to determine 
the spatial matchups of a particles with collisional depth 6S between 
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The limits of integration in (32), (35) are determined in the 
following way. The phase space map for an unscattered particle between 
"sources" r , and observer is given by (1) the pitch angle map: 

ORIGINAL PAGE Iti 
OF POOR QUALITY 

cosf0(n)) = - fj- r>) ' 


and (2) energy conservation 

l m e ) /- =* 

where 

p - ^7 


and r < refer to source regions outside (r ) or inside (r^) , that of the 
> 

observer. 

Since cos 9 must be real, the pitch angle map places spatial bounds 
on the integrals considered: 


(95 

thus implicitly defining r uS 


mt/’-fr *) £ „#■ 

Bfnj 

r L 


mv z (n, A ) _ # 

Bf%J =r 

If jB(r)j and 0 are monotonic functions there are only two roots. Energy 
conservation also places a constraint on the radial domain of integration: 


0* J- -£mv & {r*) J ^ c '$(%) 
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In practice this last constraint only impacts r^ because the particle 
is giving up kinetic energy as r increases. 

For the special case of fi* = 0, the radial limits of integration 
are determined implicitly by 

^ far*) + z e S Cr *) . - Z<r <f fa ) 

and r L is the smallest radius for which f^r^) bas sufficient population 
at kinetic energy 

jwt/ 2 fc)= is* far*) ~Zfe<jfa)) 

to be detected in comparison with the local distribution. For all the 

"realizations" which we discuss u* = Q and r T = 1.03 R . 

^ L © 


(35b) 


(35c) 
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4. IMPLEMENTATION 

In order to carry out the Integrals indicated in (35) we need to 
specify the profiles of the electron and proton temperatures, their 
density, as well as the gravitational, electric and magnetic fields. The 

4 

first three requirements define the Krook approximation to the collision 
operator through (7a) by determining f^ 0 ^ (r) and T[f^^°^ (r) ]; the last 
three profiles define the force field,’ F, of the Boltzmann evolution 
operator (1) . 

In Figure 5 we show a density profile, developed by Sittler and 

Olbert (1978), Sittler (1978) from composite experimental sources, which 

we will adopt to be the density for In Figure 6 we show the 

empirical estimates of the radial variations uf T e ^ * Kear the sun 

we have used the measurements and arguments of Bame et al. (1974), who 

—2/7 

suggest that T^ ct r ‘ near the base of the corona. We have used the 

results of Ogilvie and Scudder (1977) using data between 0.45 and 0.85 AU 

-0 3 

that T e , thermal, varies like r * and the results of Gringauz and Verigin 
(1975) for 1 < r < 1.5 AU and an asymptotic adiabatic behavior at large 
distances (> 10 AU). We have constrained the electron profile to agree 
with 1.5 x 10 J °K at 1 AU and to be l.o x 10 6 °K near 1 r q j consistent with 
coronal hole upper limits. 

To evaluate the Coulomb logarithm and assess overall energetics, a 
proton thermal profile is also required. We assume T (t) = T (r) for all 

e p 

radii inside that radial distance r = R beyond which adiabatic cooling 

4o 

determines the average < T^(l AU) > ^ 7 x 10 K. R is usually 60-90 R^, 
consistent, for example, with the two. fluid results of Acuna (1976). 

We have used the generalized Ohm's Law (c.f. e.g., Rossi and Olbert, 




T (°K ) 
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1970, p. 35Q, eqtn. (12.26)) to self-eonsistently estimate the time 
average parallel electric field in the proper frame of the plasma, E. *; 







( 36 ) 


where the additional eff ective electric field E, ( ' is given "by 




(jj?)t *■ j 



where T] is the fractional net charge density and the angular brackets 
denote a spatial average and the parenthesis ( )^ denote a temporal 
average over a multiple of fluctuating cycles. 

The last term, in (36) is non-zero when there are currents in the 
proper frame, waves propagating in the system or deviations from charge 
neutrality. In most circumstances the Hall effect term (the first one in 
(Ejj f )) is the most important. It can be shown that Alfvenic fluctuations, 
for example, give anon-zero contribution to (En’)^. Generically, (Em')^ 


represents the consequences of uncompensated forces in the proper frame 
of the plasma which have a non-vanishing time averaged component along 


< B >. 


Given the current state of our understanding of waves and their time 
averaged current systems in Che solar cavity, it is difficult to 
estimate Ey 1 (r) with precision. Nevertheless, the order of magnitude of 
the potential energy barrier which this additional electric field produces is of 
the order of the waye work necessary to have an energetically self-consistent 
solar wind expansion: 
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4 r -dh' 

AnR 0 h <r'> 

To be effective in accelerating the solar wind to the observed speeds 

this wave work is expected to be done in the transonic region of the 

inner corona- We introduce a scale, L , over which this work is 

wave 

principally done and phenomenologically have incorporated the potential 

barrier due to y W as OKIGINAL PAGE IS 

OF POOR QUALITY 


(37) 


iLe(H- + £ K)) * =4 tU/- 7 -— ‘ ’) 


/.o$£ r 


(38) 


where A ~ Work from energy conservation considerations. This portion 
of the total electrostatic barrier for electrons achieves 76% of its 
asymptotic value at = 31^^^ + 1.03 S^. Our assumption of the form 

of the additional electrostatic harrier (38) corresponds to an effective 
electric field which is strongest near the base of the corona and goes to 

-z 2 

zero as e , Z = (r -1.03 R ) /L , 

’ 0 wave 

We have therefore used an electrostatic parallel electric field of 
the form. 

v -3 a z / 

— =- + — sech ( , 

C./i, / 


= 


e d e 


^ / h- /tsfy ) 


wat/e 


(39) 


for a positive magnetic sector. 
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When solutions are reported below, I* ave » will be a parameter, but 
usually of the order of 5 1^. With these considerations,- the interplanetary 
potential is energetically self-consistent with the asymptotic wind speed 
used to define the topology of the magnetic lines of force* 

In Figure 7 we show the interplanetary potential which results from 
integrating Ey* along a typical field line with Parker topology when the 
above empirical profiles of n e (r) , T^Cr) and Ejj’(r) are used. The 
asymptotic flow speed was assumed for this example to be 4Q0 km/sec and 
ave was 5 R^. Shown for comparison is an interplanetary potential 
profile reported by Lemaire and Scherer (1971) using an exospheric solar 
wind model* 

This electrical potential represents a significant barrier for the 
escape of electrons from the proximity of the lower corona; at the same 
time it is an acceleration for the positive ions. An electric field of 
this type is implicit in fluid solar wind models which externally impose 
the local condition that. the hulk speeds and densities of the electrons 
and ions he equal (c-f., e.g., Hartle and Sturrock, 1968). This type 
electric field is also implicit in the Spitzer-Harm (1953) collision 
dominated coefficient of thermal conduction (through the assumption j* = 0) 
which is often used in solar wind models. The electric fields used in 
these references are predicated on local thermodynamic equilibrium (LTE) 
which is not justified in the solar wind. 
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5. EESDLTS 

The results of two realizations of the formal solution (15) will 
now be discussed in some detail. The two realizations differ by the 
asymptotic flow speed attained by the salat wind* but refer to the same 
spatial location of 1.1 AU. 

General Features of Solutions 

The dots in Figure 8 show a cut of the electron distribution 
function along the local magnetic field line at time, t, and position 
r = 1.1 AU, as a result of the numerical integration of (15) - (35) con- 
sidering only Coulomb collisions. The asymptotic wind speed was 400 km/sec 

and L was taken to be 5 The solid curve is the local assumed form 

wave © 

/ \ 

of, f^_ * (AS = 0, Vjj) , which is given for reference. Comparing the 

computed distribution (dots) with the solid curve and referring to (16) 
or (17) we see that the differences between the two curves are due to the 
integrated contributions originating away from the observer. 

In Figure 9 are shown fore and aft cuts along the local magnetic 
field of f^ near 1 AU reported by Feldman et al. (1975). The data of 
their Figure 14 have been transformed into the proper frame using the 
bulk velocity given in Figure 3 of that paper which was 520 km/ sec. 

Although we have not tried to fit the observed profiles by adjusting 
various parameters, we find there is a remarkable similarity between the theoreti- 
cal profile Figure 8 and the data Figure 9. We draw particular attention to the 
correspondences in (1) size, (2) shapes, and (3) skewness between the 
theoretical and experimental profiles. Of course, this is not a complete 
comparison which would require a theoretical determination of the entire 
pitch angle distribution. 



Figure 8 
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Figure 9 
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It is instructive to characterize the location where the thermal 
and extra thermal papulations at 1 AU were last members of the thermal 
collisions! population. In order to do this we recast the integrals 
of (35) as 




Lf(^\ K ?<&*) 


(40) 


-f 



Af/ e/,12, 


where Af(R^, R) is the integral over the integrands of (35) but with 
limits of R* ±AR. such that 

i i 
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In Figure 10 we have indicated the spatial location of the 
collisional distributions which dominate the populations tw0 energy 
extremes of the theoretical profile given In Figure 8. The value of 
the theoretical profile at the given energy in Figure 8 is the 
corresponding area under the indicated curves in Figure 10. 

Figure IQ clearly restates the local character of the thermal 
population. The dominant contributions for particles of 5 eV kinetic 
energy near earth come from within ±0.5 AU (approximately the thermal 
mean free path) of the observer. The forward pitch angle population 
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at 5 eV (0 = Q) come predominantly from collisional populations sunward 
of the observer; for. 0 - 180°, at 5 eV the converse is true as shown in 
the shaded contribution. 

In the extratherm al regime (as an extreme example kinetic energy 

500 eV) we see that both forward (0=0 and aft (0 = 180°) populations 

were last members of a collisional distribution deep in the base of the 

corona (< 3 and certainly near the site of the initial expansion of 

the solar wind. This is the same area where theoretical extrapolation of 

Mariner 10 radial variations of T (r) and T„(r) predicted a single 

e ti 

Maxwellian distribution. The particles observed at 500 eV propagating 

towards the sun at 1 AU (0 = 1S0°) are the hackscattered fraction of the 

outward propagating 5QQ eV particles which passed the observer, radial 

distance of the observer at some, earlier time prior to their detection. 

+ — 

The contributions to Af and Af (500 eV) are essentially in the ratio of 

6 : 1 . 

The particles with kinetic energies intermediate between these two 

extremes have collisional antecedents which vary continuously between 

the two extremes depicted in Figure 10. For some intermediate locally 

transthermal energies, there are two domains of equal importance 

(i.e., Af (R r , R) is two humped). For locally subthermal energies 

(E < IcT^) the smaller contribution is at the sun. As the local energy 

becomes increasingly extrathermal, the dominant collisional antecedent of 

the observed f recedes to lower and lower levels in the corona, 
e 

The reason for these distributed contributions to the locally 

observed f is intimately related to the spatial variation of the 
e 

interplanetary potential and the Coulomb collisional depth of the 


observer from the sun. 
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An observed electron of kinetic energy f rav (R) is energetically 
accessible from a level R' provided the kinetic energy at R- is 


i no t >/ £ f/?') -•= i m -+ Be (4 (£} - 4> fe>) ) (42) 

These accessible orbits make an observable contribution -to the electron 
distribution function at R provided that a sufficient fraction can survive 
collisional extinction en route from R 1 “+ R, as compared with the local 
collisianal Maxwellian: 


Prom Figure 7 we see that Z-(e0(R) -e0(R)) is positive (negative) 
for R 1 < R (R 1 > R) . Thus the electrons detected near earth have given 
up nearly 1 keV of kinetic energy en route from the corona to the observer. 
Electrons with energy of 500 eV at 1 AU had a kinetic energy of ~ 1.5 key 
in the coronal Maxwellian population of characteristic energy of ~ 100 eV, 
15 standard deviations off of the maximum phase density at R* ~ 1-2 R^. 
because the relative density of the lower corona becomes large so rapidly 
and because 6s (R 1 , R, 500 eV) is so small Af(R f , R) can give a measurable 
phase density signal to the observer near 1 AU. 

On the other hand because the Coulomb collisional extinction 
6S, depends so strongly on the kinetic energy of the particle (c.f., 

(10a, 12)) it is very difficult to get non-local contributions to f 
observed for subthermal energies (E <kT c (R)). Although condition 

(42) is readily met at a reasonable range of non— local R 1 s, the survival 
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probabilities are so small (SS (R , R) for kinetic E<k^T are so large) 

that only those collisional distributions within the vicinity can survive 

to the observer. [Mathematically, 6S (R* , R) -> 0 for any E provided 

R’ (t T ) -* R(t) is small enough: c.f., equation (12).] 

We thus see that the population of the forward extrathermals should 

in steady state be a very sensitive indicator of the interplanetary 

potential barrier between the base of the corona and the observer. The 

best local indicator of a large potential barrier to the corona for 

electrons is a higher solar wind speed, other things being equal such as 

the V*Pe contribution to the potential drop. 
nTr) 

The above conclusion implies an anti-correlation between local plasma 
bulk speed (if its nearly the asymptotic bulk speed) and the local 
fractional extrathermal (n^^/n) phase density. From our discussion of 
Figure 10 we should not interpret this as a local cause and effect 
relationship, rather — the cause of both local effects is in the strong 
acceleration region removed by 1 AU from the observer. This global sense 
of causality in the inner heliosphere must prompt reexamination of the 
use of transport coefficients and equations of state which by definition 
relate local forces to local responses. We return to these matters in 
detail in paper III. 

In Figure 11 we show in the same format as in Figure 8 a solution 

for U_ ■= 800 km/sec with L still = 5 R . and all other parameters 
03 wave © 

which are independently specifiable the same. As explained previously, 
the net wave work is adjusted to give consistency to the expansion energetics 
and magnetic topology. As can he readily seen the suprathermal population 
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is depleted in a relative sense as we qualitatively expected. Further- 
more, the estimate of the "differential temperature" for the suprathermals 
is also reduced when the asymptotic flow speed is increased. The heat 
flux {q^l depends on the number of carriers available and it is known that 
the most important contribution to [q [ at 1 AU occur between 25 eV and 
130 eV (Scudder, 1970 and Montgomery, 1970) centered on the break region. 

If we thus deplete this population and lower its most probable speed 

h 

(c&T the heat flux it carries must also be reduced. We thus expect 
on the basis of this model that the interplanetary heat flux density 
in quasi-steady state flows should be anti-correlated with the local 
bulk speed. As before this is a global cause with two local signatures; 

[u{ and [q | . There is thus no intrinsic local thermodynamic inter- 
relationship between Jq g | and J U | ; rather, there is a global integro- 

differential interrelationships hetween the whole expansion history of the flow 
tube and all pressure and density gradients (which determine Ej|) along it. 

Figure 12, shows the dominant spatial contributions to the 
5 and 500 eV points of Figure 11- Note the changed scale of the left 
ordinate relative to Figure (10). In this higher speed Cu^ = 800 km/sec) 
example, the deep corona is essentially inaccessible to the normal halo 
population. The sparser contributions come from 8-20 with somewhat cooler 
collisional populations which are sufficiently populated at the correct 
shifted energies to have energetic access and phase density excess 
sufficient to overcome extinction. The contributions to the local population 
of f remain essentially immune to the change in the global dynamics, 
which in this example merely changes the floating potential of the 
vicinity of the observer hy a constant amount. 
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Thermal Domain Effects of Scattering 


Returning to Figures 8 and 11 and concentrating on the lower energy 
portion of the distribution we see on first impression that the results 
of the scattering have not changed the local distribution. However, 
the slight departures of the dots from solid AS = 0 profiles in Figures 8 and 11 


are significant and reflect the fact that so long as the mean 
free path of the medium is finite there will be corrections to the 
asymmetries of the resulting f (v) . These low energy modifications are 
most nearly in response to local thermodynamic forces as will be discussed 
in detail in 1X1. 

Solar Rotation Effects 

The asymptotically wrapped archimedian spiral of the magnetic tubes 
of force in the plane of the ecliptic plays an important role in de- 
termining the nature of f in the ext rath ermal energy regime reported in 

the previous solutions. The arc length JL along B scales (c.f. , Equation C.345) 
2 

as r for large r in the ecliptic, but only linearly with r oyer the 
magnetic poles of the sun. Thus the cumulative number of Coulomb '’collisions" 
is relatively enhanced in the ecliptic versus the polar regions by this 
geometrical hindrance of the guiding centers following the tubes of force. 

This situation may have an important impact on the nature of the magnitude 
of Coulomb back scattered extrathermals seen over the pole of the sun. 

Whether this implies a large skewness and resulting large heat flux on polar 
tubes of force depends on the electrostatic shielding of the interplanetary 
medium from the corona. The magnitude of this shielding is determined by 
the detailed polar coronal profiles of density, electron pressure and 
waves. If the polarization potentials are very strong so as to support large 


42 


asymptotic flows over the poles, there may not he a substantial extra- 
thermal population "leaking” into interplanetary space. In this situation 
it may not he important for the value of the heat flux that cumulative 
Coulomb hackscattering of extrathermals is weaker than In the ecliptic. 
Completion of the Iterative Scheme of Equation (7) 

Cl) 

The solutions, f , presented in this section represent the first 
step of the chain of iterations outlined in Equation (7) . These first 
iterates already very closely resemble typical data. From those equations 
and Equation (35) we see, however, that to perform the next iteration to 

obtain we must first determine £^(r, v) at all positions along the 

( 2 ) 

tube of force since f is an integral over all screened phase densities 
that have access to the observer. In addition, the collision time t is 
also a functional of f^(j:, v) . Since f^ (r. v) Is only available 
numerically the velocity dependence of T would also have to he determined 
numerically. This is not attractive (c.f., Rossi and Olbert, 1970, 

Equation 11.40-11.45). Subsequent stages of the iteration process will 
become more and more intractable. 

The solutions for are not the fully self-consistent solution to 

the Coulomb scattering problem in the solar wind. These solutions indicate 
what type of modification the finite Coulomb mean free paths have on 
shaping the local distribution function and what forces and scales are 
important in the problem. We shall formulate and solve the more difficult 
fully self-consistent problem in paper IV of this series. 
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6. conclusions ORIGINAL PAGE IS 
OP POOR QUALITY 

We have shown that Coulomb collisions are substantial mediators 
of the interplanetary electron velocity distribution function. This 
is true both for the thermal distribution and the suprathermal 
populations. For the suprathermals, the cumulative effects of Coulomb 
interactions take place on the scale of the heliosphere itself, whereas 
the thermal population’s Coulomb interactions are numerous on the local 
scale (% - % AU) near the point of observation (1 AU). The suprathermal 
electron phase density at 1 AU is populated by electrons which most 
recently were members of collisional distributions very deep in the 
dofona in or near the region (1-10 of the strong acceleration of the 
solar wind. 

We have also shown that the properties of Coulomb collisions and 

the forces implied in the solar wind expansion, without benefit of 

wave-particle effects, can place a zone for a bifurcation of the electron 

distribution function deep in the corona. The site of this fractionation 

(1-10 R q ) is consistent with the result of extrapolating the radial profiles 

of T and T„ determined from Mariner 10 data reported by Ogilvie and 
c H 

Scudder (1978) . 

In the steady state the present considerations lead to a variety of 
predictions which can be tested against the growing body of electron 
data becoming available. Some of the more obvious relations implied 
in the model for an observer of a fixed radial position are the anti- 
correlation of local bulk speed with fractional extrathennal density, 
extrathermal differential temperature estimates, and heat flux density. 

In addition the extrathermal differential temperature should he nearly 
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independent of heliocentric distance within 1 AD due to their small 
number of Coulomb momentum transfer ’’collisions” and thus even fewer 
energy exchange ’’collisions” between the corona and 1 AU. 

We strongly emphasize the need to think of the trans and extra- 
thermal electron kinetics in the solar wind on the global scale: 
because 1) the causes of in situ changes in properties of these electrons 
may not always be local ones; 2) non-local causes may affect different 
local parameters (e.g., speed, suprathermal population) in different 
ways, thereby inducing correlations in local measurements and the con- 
verse: 3) observed local correlations in electron parameters do not 
imply local causal interrelationships; and lastly 4 } the conventional 
transport question for solar wind electrons, must be reposed to reflect 
the global character of the trans and suprathermal electrons which 
implement the various transport signatures. 

By contrast, th i low-energy electrons should in this picture 

reflect local cause and effect relationships. These electrons should 

be shaped by the local thermodynamic forces. Nevertheless, the thermals and 

sup rathe rmals are Coulomb collisionally interacting locally- In steady state, the 

thermal and suprathermal electrons observed at the same time were once members 

at different times of the (same) collisional equilibrium at the base of 

the same field line. Thus in varying degrees, the halo more, the core 

less, the sub and suprathermals retain a memory of having interacted 

with this same reservoir. We thus suggest that the nearly fixed ratio 

between T and T _ at 1 AU, for example, is a remnant of this fact, 
e extra 

This is another example of a non-local factor inducing local correlations 
in observed parameters. 
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ORIGINAL PAGE IS 
OP POOR QUALITY 

By considering the global effects of Coulomb collisions, we have 
shown that the distribution function for electrons simij ar to those 
reported can be reproduced in overall shape, skewness, and partition of 
phase density with energy. In particular we have no difficulty in 
obtaining the backscattered portion of the extrathermal population 
(the "halo"). The extrathermal population is a natural consequence of 
the evolution of a very hot, dense coronal plasma on open field linos. 

To summarize, we should like to re-emphasize the point made in 
the introduction: the observed "overabundance" of higher-energy 

(sup rathe rmal) electrons as judged by the local Maxwellian distribution 
of lower-energy (core) electrons are a direct con- 
sequence of the non-uniformity of the expanding medium of wftu:h the 
electrons are an intrinsic part. Expansion implies cooling and decrease 
in density as one moves along any open magnetic tube of force. This 
circumstance causes the Coulomb mean-free-path of sufficiently high- 
energy electron to grow with distance to he point that it becomes 
comparable or larger than the local scale length; this, in turn, renders 
the problem to be global (non-local) in character. The growth o £ the 
mean-free path does not, however, imply that the effects of Coulomb 
collisions become negligible on a global scale. In fact, for the 
topology of magnetic fields at hand, the Coulomb collisions are never 
negligible, whatever the energy of the observed electron may be. 

The electrostatic magnetic and gravitation fields combined with the 
action of Coulomb processes on global scale shape a 
non-Maxwell ian distribution that, at first glance, appears to be 
composed of two distinct electron populations, thermal and suprathermal. 
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Unlike the customary discussion of the origin of suprathermals (Parker 
.and Tidman, 1958). We do not invoke Fermi, or betatron equivalent 
energisation nor have we invoked shock acceleration, nor do we require 
leakage from closed coronal field lines as suggested by Feldman eit al. 

(1975) to populate the extrathermal regime. 

In conclusion, we would like to stress that the above described 
set of circumstances is not peculiar to our sun alone, but rather obtains in 
a great variety of ^ro tat ingj astro physical objects. Any star (or, even 
perhaps some planets like Jupiter) that possesses an expanding ionized 
atmosphere with open magnetic field line topology should have electron 
distributions of the type observed in our solar system. 

The above discussion should be not misconstrued as to indicate that 
the interactions of positive ions with waves are ruled out by this 
analysis. In fact, they may be of primary importance. Provided 
the magnetic turbulence invoked for the ions is laminar on the scale 
of the electron gyroradius there will be no contradictions to the 
current work since such variations would appear to make only minor path 
length modifications to 6S (R 1 , R) , without there being any significant 
electron -wave scattering interaction. 

We believe that our results are sufficiently encouraging to warrant 
a more detailed and self-consistent formulation of the Coulomb scattering 
problem in the solar wind which will be discussed in paper IV of this 
series. 
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FIGURE CAPTIONS 


Figure 1 Idealized magnetic topology. For electrons the arc 

length, time coordinates (if.*, t*) are synonynous with 
the space time position (r*, t*) , 


Figure 2 


(a) The first five types of space-time trajectories with 


* 


access to (&*, t*) with v^ 0; (b) first five types of 

space-time trajectories with access to the same (&*, t*) 

A 

position but with v,. = - v«. 


Figure 3 Schematic of fore, f , and aft, f , electron distribution 

functions under different assumptions of the dominant 

manner in which suprathezmals are hackscattered. A 

Maxwellian in this format is a straight line. (A) Exclusively 

Elastic: in this limit f /f -r for E > 7 kT and 

o c 

— d d f in extrathermal regime. (B) Elastic 

Dominant plus Inelastic ; (-^/d bn f*) > -(^/d bn f ) 
extrathermal; (-^d to f") extra > & dermal* 


(C) Inelastic Exclusive: (-^^d bn f + ) > (-^*/d to f ) 


extra 


but ) ext ra~ 

kT(R*) where R* is the distant radial position where extra- 
thermals have lost so much kinetic energy as to become a 
member of the local Maxwellian with temperature ki(R*) • 




Figure 4 The ratio of fore to aft phase density (f /f ) as a function 

of proper frame velocity parallel to B.as adapted from published 

data open circles Feldman _et al . , 1975 and solid dots from 

illustrative theory— ‘this paper Figure 8. Horizontal line gives 

theoretical asymptote of which prevails for excusively 
elastic case of Figure 3a„ 

ORIGINAL PAGE IS 
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Figure 5 


Figure 6 


Figure 7 


Figure 8 


Figure 9 


Figure 10 


Radial density profile after Sittler and Olbert (1978) 
and Sittler (1978). 

Radial thermal electron and proton temperature profile. 
Constructed as described in text based on empirical 
constraints. 

Self-consistent change in potential energy for an electron 
where the zero of potential has been shifted to the base 
of the corona. An electron gives up kinetic energy which 
is stored as this potential change in leaving the lower 
corona. 

, . ( 1 ) * * 

(a) Theoretical profile of f v * (vy . R ~ 1.1 AU) for 

U OT = 400 km; (solid curve) f (AS = 0) used in solution; 

see text for details. 

Observed f (vy , r ^ 1 AU) adapted (unfolded) from 
Feldman et_ al_. (1975) Figure (1 ). 

Spatial locations (R T ) of the dominant contributions 

* * A (1) * 

Af(R , R j vy , E ) which determine f (vy , R jin Figure 8. 

The low energy (5 eV) thermal population were last members 

of a collisional distribution very near the observer. The 

high energy (500 eV at 1 AU) supra thermals (halo) 9=0 and 

0 = 180° were last members of a collisional distribution 

deep in the corona. The 0 = 180° halo particles are 


scattered backwards toward the observer at R' ~ 2-10 AU 
which is outside the observer's position of ~ 1.1 AU. 



Figure 11 


Figure 12 


Theoretical profile (?) (vy , R ^ 1.1 ATT) for 

= 800 km/sec. Solid curve is f^ 0 ^ (vy , AS - 0); 
see text. 

Dominant spatial contributions to solution in Figure 11, 
same format as Figure 10, but for the solution shown in 
Figure 11. Kote reduced scale for extrathermal contri- 
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